Cosmological dynamics of dark matter Bose-Einstein Condensation 
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Once the critical temperature of a cosmological boson gas is less than the critical temperature, a 
Bose-Einstein Condensation process can always take place during the cosmic history of the universe. 
In the Bose-Einstein Condensation model, dark matter can be described as a non-relativistic, New- 
tonian gravitational condensate, whose density and pressure are related by a barotropic equation of 
state, with barotropic index equal to one. In the present work, we study the Bose-Einstein Conden- 
sation process in a cosmological context, by assuming that this process can be described (at least 
approximately) as a first order phase transition. We analyze the evolution of the physical quantities 
relevant for the physical description of the early universe, namely, the energy density, temperature 
and scale factor, before, during and after the Bose-Einstein Condensation (phase transition). We 
also consider in detail the epoch when the universe evolved through a mixed condensate - normal 
dark matter phase, with a monotonically growing Bose-Einstein dark matter component. An im- 
portant parameter characterizing the Bose-Einstein Condensation is the condensate dark matter 
fraction, whose time evolution describes the time dynamics of the conversion process. The behav- 
ior of this parameter during the cosmological condensation process is also analyzed in detail. To 
study the cosmological dynamics and evolution we use both analytical and numerical methods. The 
presence of the condensate dark matter and of the Bose-Einstein phase transition could have mod- 
ified drastically the cosmological evolution of the early universe, as well as the large scale structure 
formation process. 
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I. INTRODUCTION 

At very low temperatures, particles in a dilute Bosc 
gas can occupy the same quantum ground state, forming 
a Bose-Einstein (BEC) condensate, which appears as a 
sharp peak over a broader distribution in both coordi- 
nates and momentum space. The possibility to obtain 
quantum degenerate gases by a combination of laser and 
evaporative cooling techniques has opened several new 
lines of research, at the border of atomic, statistical and 
condensed matter physics (for recent reviews on the Bose- 
Einstein Condensation see [H-Q)- 

An ideal system for the experimental observation of the 
Bose-Einstein Condensation is a dilute atomic Bose gas 
confined in a trap and cooled to very low temperatures. 
BEC were first observed in 1995 in dilute alkali gases 
such as vapors of rubidium and sodium In these ex- 
periments, atoms were confined in magnetic traps, evap- 
oratively cooled down to a fraction of a microkelvin, left 
to expand by switching off the magnetic trap, and sub- 
sequently imaged with optical methods. A sharp peak 
in the velocity distribution was observed below a critical 
temperature, indicating that condensation has occurred, 
with the alkali atoms condensed in the same ground state. 
Under the typical confining conditions of experimental 
settings, BEC's are inhomogeneous, and hence conden- 
sates arise as a narrow peak not only in the momentum 
space but also in the coordinate space 
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Since in the terrestrial experiments Bose-Einstein Con- 
densation is a well-known phenomenon, the possibility 
that a similar condensation may have occurred during 
the cosmological evolution of the universe cannot be ex- 
cluded a priori. In fact, once the critical temperature 
of the boson gas is less than the critical temperature, 
BEC can always take place at some moment during the 
cosmic history of the universe. Different aspects of the 
BEC cosmological transition were considered in 
The critical temperature for the condensation to take 
place is T cr < 2irh 2 n 2 / 3 /rnks, where n is the particle 
number density, m is the particle mass, and ks is Boltz- 
mann's constant [J-Q- On the other hand, cosmic evo- 
lution has the same temperature dependence, since in 
an adiabatic expansion process the density of a matter 
dominated universe evolves as p oc T 3 / 2 171. Therefore, 
if the boson temperature is equal, for example, to the 
radiation temperature at a redshift z = 1000, the criti- 
cal temperature for the Bose-Einstein Condensation is at 
present T cr = 0.0027-ftT Q. Since the matter tempera- 
ture T m varies as T m oc a -2 , where a is the scale factor 
of the universe, it follows that during an adiabatic evo- 
lution the ratio of the photon temperature T 7 and of the 
matter temperature evolves as T 7 /T m oc a. Using for 
the present day energy density of the universe the value 
p cr = 9.44 x 10 -30 g/cm 3 , BEC takes place provided that 
the boson mass satisfies the restriction m < 1.87 eV [f|. 
On the other hand, we expect that the universe is always 
under critical temperature, if it is at the present time [7J . 

Despite the important achievements of the A Cold 
Dark Matter (ACDM) model @, at galactic scales ~ 10 
kpc, the standard cosmological model meets with severe 
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difficulties in explaining the observed distribution of the 
invisible matter around the luminous one. In fact, TV- 
body simulations, performed in this scenario, predict that 
bound halos surrounding galaxies must have very charac- 
teristic density profiles, called the Navarro-Frenk- White 
(NFW) profiles, which feature a well pronounced cen- 
tral cusp PNFw{ r ) = Ps/(r/r s )(l + r/r s ) 2 , where r s is a 
scale radius and p s is a characteristic density [1] . On the 
observational side, high-resolution rotation curves show, 
instead, that the actual distribution of dark matter is 
much shallower than the above, and it presents a nearly 
constant density core (Io| . The core-cusp problem of 
standard CDM models can be solved by assuming that 
dark matter is composed of ultralight scalar particles, 
with masses of the order of 10~ 22 eV, initially in a cold 
Bose-Einstein condensate ( "fuzzy dark matter") . 
The wave properties of the dark matter stabilize gravita- 
tional collapse, providing halo cores. On the other hand 
dark matter models with pressure, satisfying a polytropic 
equation of state, give an excellent fit to the observed 
galactic rotational curves [l^ . The polytropic equation 
of state can describe extended theories of dark matter in- 
volving self-interaction, non-extensive thermostatistics or 
boson condensation (in a classical limit). In such mod- 
els, the flat-cored mass profiles widely observed in disc 
galaxies are due to innate dark physics, regardless of any 
baryonic motion. The fine-structure in the observed inner 
mass distribution of the Milky Way can be explained only 
if the infalling dark matter particles, from which such 
galactic halos formed, had a net overall rotation, causing 
a tricusp caustic ring of dark matter [l3[ . For standard, 
non-interacting CDM models, however, one expects the 
infall to be irrotational. Bose-Einstein condensate dark 
matter can form vortices, thus leading to a net overall 
rotation of the galactic halos. Hence Milky- Way observa- 
tions may have already detected some specific signatures 
of the Bose-Einstein Condensate dark matter [1J| . 

The possibility that the galactic dark matter is in the 
form of a (cold) Bose-Einstein condensate was considered 
in detail in [15]. The density distribution p of the static 
gravitationally bounded single component dark matter 
Bose-Einstein condensate is given by p (r) = p c sin kr/kr, 
where p c is the central density of the condensate, p c — 
p(0), and k is a constant. At the boundary of the 
dark matter distribution p(R) = 0, giving the condition 
kR = 7r, which fixes the radius of the condensate dark 

matter halo as R = tt^J h 2 l a /Gm 3 , where l a is the s-wave 
scattering length, and m is the particle mass. The total 
mass of the condensate dark matter halo M can be ob- 
tained as M = 47r 2 (h 2 l a /Gm 3 ) 3/2 p c = \R?p c j-K, giving 
for the mean value < p > of the condensate density the 
expression < p >= 3p c /7r 2 . The mass of the particle in 
the condensate is given by 

rn = j « 6.73 x 1(T 2 x 

[Ufm)] 1/3 [i?(kpc)r 2/3 eV. (1) 



The global cosmological evolution and the evolution of 
the density contrast in the Bose-Einstein condensate dark 
matter model, in the framework of a post-Newtonian cos- 
mological approach, were investigated in [l6[ and [l7| . 
The global cosmological evolution as well as the evolu- 
tion of the perturbations of the condensate dark matter 
shows significant differences with respect to the pressure- 
less dark matter model, considered in the framework of 
standard cosmology. Therefore the presence of conden- 
sate dark matter could have modified drastically the cos- 
mological evolution of the early universe, as well as the 
large scale structure formation process. The core/cusp 
problem and the dark halo properties of the dwarf galax- 
ies were recently analyzed in [18| . Different properties of 
the dark matter condensate have been extensively stud- 
ied in the physical literature 



It is the purpose of the present paper to investigate 
the cosmological Bose-Einstein Condensation process of 
dark matter. The condensation process is interpreted as 
a phase transition taking place sometimes during the cos- 
mic history of the universe. A Bose-Einstein Condensa- 
tion process of the dark matter in the expanding universe 
can be described generically as follows. As the normal 
bosonic dark matter cools below the critical temperature 
T cr of the condensation, it becomes energetically favor- 
able to form a condensate, in which all particles are in 
the same quantum state. However, the new phase does 
not appear instantaneously, and the two phases coexist 
for some time. The transition ends when all normal dark 
matter has been converted to a condensed state. During 
the transition phase the global cosmological evolution of 
the universe is changed. By assuming that the phase 
transition is of the first order, we study in detail the evo- 
lution of the relevant cosmological parameters (energy 
density, temperature, scale factor, etc) of the normal dark 
matter and Bose-Einstein condensed dark matter phase, 
and the condensation process itself. An important pa- 
rameter to describe the Bose-Einstein Condensation is 
the condensate dark matter fraction, whose time evolu- 
tion describes the time dynamics of the conversion pro- 
cess. The behavior of this parameter is also analyzed in 
detail. 



The present paper is organized as follows. In Sec- 
tion [TT1 we briefly outline, for self-completeness and self- 
consistency, the basic properties of the normal and Bose- 
Einstein condensed dark matter. We also lay down the 
equations of state and the relevant physical quantities 
that are analyzed in the remaining Sections. The cos- 
mological dynamics of the Bose-Einstein Condensation 
of dark matter is analyzed in Section Hill We discuss and 
conclude our results in Section IIV1 In the present paper 
we use the CGS system of units. 
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ds 2 = -c 2 dt 2 + a 2 (t) (dx 2 + dy 2 + dz 2 ) 



II. COSMOLOGICAL DARK MATTER IN 
NORMAL AND BOSE- EINSTEIN CONDENSED 
STATES 

In the present Section we outline the relevant physical 
processes and quantities of the Bosc-Einstein Conden- 
sation, which will be used in the following Sections to 
study the dynamics of the condensation and some of its 
cosmological implications. 

We assume that the space-time geometry is the fiat 
Friedmann-Robertson-walker (FRW) metric, given by 

(2) 

where a is the scale factor describing the cosmological 
expansion. For the matter energy-momentum tensor on 
the we restrict our analysis to the case of the perfect fluid 
energy-momentum tensor, 

T^ v = (pc 2 + p)u^u v + pgH". (3) 

The Hubble function H(a) is defined as H — a/a. As 
for the matter content of the universe, we assume that it 
consists of radiation, with energy density p ra d and pres- 
sure Prad, pressureless (pb = 0) baryonic matter, with 
energy density pb, and dark matter, with energy den- 
sity p x , and pressure p x , respectively. In the following 
we neglect any possible interaction between these com- 
ponents, by assuming that the energy of each component 
is individually conserved. Thus, the gravitational field 
equations, corresponding to the line element @ become 

A 2 



8TTG (p b + p r ad + p x ) + A, 



8ttG 



Pi 



{Pb +Prad +P X ) + ^ 



(4) 



(5) 



Pi + 3 ( pi + ^ ) - = 0, i = b, rad, \- (6) 
\ cr / a 

The cosmological evolution of the energy density of 
the baryonic matter and radiation are given by pb = 
pbfi/ (a/ao) and p ra d = p ra d,o/ (a/a ) , respectively, 
where pb,o and p ra d,o are the energy densities of the mat- 
ter corresponding to the value a = ao of the scale fac- 
tor. For the dark matter we consider a general density 
evolution of the form p x = p x ,o/ f (a/ao), where p Xt o is 
the value of the energy density of the dark matter at 
a = ao, and / (a/ao) is an arbitrary function of the 
scale factor, depending on the particular dark matter 
model. By introducing the critical density at a = ao 
as p cr fl = 3H 2 /8ttG, where H = H (ao), and the den- 
sity parameters fi^o = Pi,o/Pcr,o, i — b,rad,\, we obtain 
the basic equation describing the dynamics of the cosmo- 
logical models as 



a - H-2 
a z 



n 



6,0 



^rad,0 



n 



(a/ao) 3 (a/ao) / ( a /°o) 



(7) 



where £Ia is the density parameter of the dark energy. 
The density parameters satisfy the relation fib o+^rad o+ 



A. Normal dark matter in the early universe 

We assume that in the early stages of the evolution of 
the universe dark matter consisted of bosonic particles of 
mass m x and temperature T, originating in equilibrium 
and decoupling at a temperature Tn or chemical poten- 
tial p » m x . By assuming that the dark matter forms 
an isotropic gas of particles in kinetic equilibrium, the 
spatial number density is given by 



J 4nf(p)p 2 dp, 



(8) 



where h is Planck's constant, g is the number of helicity 
states, and 



f(p) = {exp [(E - n)] - 1} \ 



(9) 



where p is the momentum of the particle and E = 
^Jp 2 + m x c 4 is the energy. A particle species that de- 
couples in the early universe from the remaining plasma 
at temperature To redshifts its momenta according to 
p(t) = pDao/a(t), where a(t) is the cosmological scale 
factor and ao is the value of the scale factor at the de- 
coupling. The number density of the particles n evolves 
as n x ~ a~ 3 (t) [20(. The distribution function / at a time 
t after the decoupling is related to the value of the distri- 
bution function at the decoupling by f(p) — f (pa/ao)- 
The distribution function keeps an equilibrium shape 
in two regimes. In the extreme-relativistic case, when 
E pa pc, T = Tnan/a, and p = pdo>d/o>, respec- 
tively, the distribution function is given by Jer(p) — 
{exp [(pc — p)] — 1} .In the non-relativistic decoupling 
case E — p rj p 2 /2m x — pkin, where we have defined 
Pkin = P — Tix c2 ' ^ ne distribution function is Jnr(p) = 
{exp [(p 2 /2m x — Pkin)] — l} • In the non-relativistic 
case pkin and T evolve as pun = Vkin,D (an/a) 2 and 
T = T D / (a D /af, respectively 

The kinetic energy-momentum tensor associated to 
the frozen distribution of dark matter is given by 



9 / P^Pu 



^f(p)d A p. 



v h 3 J p° 
The energy density e of the system is defined as 



' = 3^8 / Ef(p)d 3 p, 



(10) 



(11) 



n 



Q A = 1. 



while the pressure of a system with an isotropic distribu- 
tion of momenta is given by 

r = i*Jpvm<i 3 p = isJ C ^ff(p)^ (12) 

where the velocity v is related to the momentum by v = 
pc 2 / E [21]. In the non-relativistic regime, when E rj 
m x c 2 and p ~ m x v x , the density p x of the dark matter 
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is given by 
obtained as [2C 



m x n x 



while its pressure P x can be 



J < i3) 



giving 



P x = p x c 2 a 2 , 



(14) 



where a 2 — (v 2 )/3c 2 , and (v 2 ) is the average squared 
velocity of the particle, a is the one-dimensional veloc- 
ity dispersion. The cosmological dynamics of the dark 
matter density is described by the equation 



p x + 3p x (l + a 2 )- = 0, 

with the general solution given by 

Px,o 



Px 



(a/ao 



(15) 



(16) 



where p x o is the density of the dark matter at a = clq. 
The dynamics of the universe with normal dark matter 
component is described by the equation 



= m 



n 



?>.o 



rad.O 



n 



xfi 



(a/aof (a/a ) 4 (a/a ) 3(1+CT2) 

In the following we consider that a = ao is the present 
day scale factor. Therefore for the Hubble constant we 
adopt the value H = 70 km/s/Mpc = 2.27 x lCT 18 s -1 , 
giving for the critical density a value of p cr ,o = 9.24 x 
10~ 30 g/cm 3 . The present day dark matter density pa- 
rameter is S7 Xi0 fa 0.228, Qb,o = 0.045, fl ra d,o = 8.24 x 
10~ 5 , and Qa = 0.73, respectively [22]. We also introduce 
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S. 



the Hubble time, defined as tn = 1 /Ho — 4.39 x 10 
Since dark matter is non-relativistic, the global cosmo- 
logical evolution of the universe is very little influenced 
by the variation of the numerical values of a 2 . 



B. Bose-Einstein condensed dark matter 

At very low temperatures, all particles in a dilute Bosc 
gas condense to the same quantum ground state, form- 
ing a Bose-Einstein Condensate (BEC). Particles become 
correlated with each other when their wavelengths over- 
lap, that is, the thermal wavelength At is greater than 
the mean inter-particles distance I. This happens at a 
temperature T cr w 2tt x h 2 p 2 l z /m 5 ^ 3 kB, where m is the 
mass of the particle in the condensate, p is the density, 
and ks is Boltzmann's constant [l|. A coherent state de- 
velops when the particle density is enough high, or the 
temperature is sufficiently low. We assume that the dark 
matter halos are composed of a strongly - coupled dilute 
Bose-Einstein condensate at absolute zero. Hence almost 
all the dark matter particles are in the condensate. In 



a dilute and cold gas, only binary collisions at low en- 
ergy are relevant, and these collisions are characterized 
by a single parameter, the s-wave scattering length l a , 
independently of the details of the two-body potential. 
Therefore, one can replace the interaction potential with 
an effective interaction Vj (r ' — r) = X5 (r ' — r), where 
the coupling constant A is related to the scattering length 
l a through A = Airfi 2 l a /m x [l|. The ground state prop- 
erties of the dark matter are described by the mean-field 
Gross-Pitaevskii (GP) equation. The GP equation for 
the dark matter halos can be derived from the GP en- 
ergy functional, 



Em = J 



'Wr^ + ^hMrOI 4 



2m 



dr 



dfdr', (18) 



where ip (r) is the wave function of the condensate, and 
Uo = 4:Trh 2 l a /m x [l|. The first term in the energy func- 
tional is the quantum pressure, the second is the interac- 
tion energy, and the third is the gravitational potential 
energy. The mass density of the condensate dark matter 
is defined as 

p x (r) = m x \if) (r)| 2 = m x p(r,t) , (19) 

and the normalization condition is N = J \tp (r)\ 2 dr, 
where N is the total number of dark matter particles. 
The variational procedure SE [ip] — fiS J \ip (r)\ dr = 
gives the GP equation as 



2m 



-S7 2 ip{r) +m x V (r)ip(r) + 



U \ip (r)\ i> (r) = pip {f) 



(20) 



where p is the chemical potential, and the gravita- 
tional potential V satisfies the Poisson equation \7 2 V — 
AirGp. In the time-dependent case the generalized Gross- 
Pitaevskii equation describing a gravitationally trapped 
rotating Bose-Einstein condensate is given by 



. d ._ . 
in—tp{r,t) = 



2m, 



■m x V(r) + U \^(r,t)f 



x 

iP(r,t). 



(21) 



The physical properties of a Bose-Einstein condensate 
described by the generalized Gross-Pitaevskii equation 
given by Eq. (|2ip can be understood much easily by us- 
ing the so-called Madelung representation of the wave 
function [1-3], which consist in writing ip in the form 



V> (r, t) = yj p (f, t) exp 



- n sm 



(22) 



where the function S (f, t) has the dimensions of an ac- 
tion. By substituting the above expression of ip (r, t) into 
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Eq. (12TT) . it decouples into a system of two differential 
equations for the real functions p x and v, given by 



(23) 



- + (,.V), 



V 



Wq, (24) 



where we have introduced the quantum potential Vq = 
— (h 2 /2m x ) V 2 A /Px/ V / Px, and the velocity of the quan- 
tum fluid v = VS/m x , respectively. The effective pres- 



sure of the condensate is given by 



Pv 



Px 



^oP x , 



where 



u 



2irTi l a 



(25) 



(26) 



The Bosc-Einstein gravitational condensate can be de- 
scribed as a gas whose density and pressure are related 
by a polytropic equation of state, with index n = 1 
15]. When the number of particles in the gravitation- 
ally bounded Bose-Einstein condensate becomes large 
enough, the quantum pressure term makes a significant 
contribution only near the boundary of the condensate. 
Thus the quantum stress term in the equation of motion 
of the condensate can be neglected. This is the Thomas- 
Fermi approximation, which has been extensively used 
for the study of the Bose-Einstein condensates [lrl3| ■ As 
the number of particles in the condensate becomes in- 
finite, the Thomas-Fermi approximation becomes exact. 
This approximation also corresponds to the classical limit 
of the theory. From its definition it follows that the veloc- 
ity field is irrotational, satisfying the condition V x v = 0. 

The cosmological evolution of the energy density of the 
Bose-Einstein condensate is determined by the equation 



Px 



30* (1 + 3 



with the general solution given by 



Px = 



(a/a ) -(u /c 2 )C x 



(27) 



(28) 



where C x is an arbitrary constant of integration. By 
using the condition that p x — p Xt o for a = do, we obtain 
the density of the condensate in the form 



Pox 



Px = 

wo {a/a o y - p 0x 

where we have denoted 



Pox 



P x ,qu /c 2 _ pcr,a^ X fl u o/ c2 

1 + p xfi U /c? 1 + /5 C r,oS\,oWc 2 



(29) 



(30) 



III. COSMOLOGICAL DYNAMICS OF THE 
BOSE-EINSTEIN CONDENSATION 

The order of the phase transition of the Bose-Einstein 
Condensation in an interacting Bose systems, i.e, the pas- 
sage from the normal to the condensed phase where all 
particles occupy a single-particle state, has been inten- 
sively discussed in the recent physical literature. Ac- 
cording to some results the Bose-Einstein Condensation 
shows a spontaneous U(l) gauge symmetry breaking, 
with the condensate fraction Nq/N playing the role of the 
order parameter [23[ . That would suggest a second order 
phase transition. However, by analyzing the dependence 
of the chemical potential p on the temperature T and 
particle density p in the framework of several theoretical 
models, describing the thermodynamical transition from 
the normal to the BEC phase in weakly interacting Bose 
gases, it was shown that none of them predicts a second- 
order phase transition, as required by symmetry-breaking 
general considerations [23 ]. These results would imply 
that the Bose-Einstein Condensation represents a first 
order phase transition 0, On the other hand, by 
definition, first-order phase transitions in thermodynam- 
ics feature a genuine mathematical singularity. Whether 
finite systems in nature can literally exhibit such a be- 
havior is a long-standing controversial question in physics 
[26| . However, a rigorous investigation into the thermo- 
dynamic instability of an ideal Bose gas, confined in a 
cubic box has shown that a system consisting of a finite 
number of particles can exhibit a discontinuous phase 
transition that features a genuine mathematical singu- 
larity, provided that not the volume, but the pressure 
is kept constant (25[. This result was obtained without 
assuming a thermodynamic limit or a continuous approx- 
imation. Hence, it seems presently that the best possible 
description of the intrinsic dynamics of the Bose-Einstein 
Condensation can be obtained within the framework of 
a first order phase transition. 



A. Cosmological parameters at the condensation 
point 

Generally, the chemical potential p of a physical sys- 
tem is a function of the temperature and particle density 
n = N/V, p = p(n,T), where N is the total particle 
number, and V is the volume. The extensivity prop- 
erty of the Helmholtz free energy F = F(N,V,T) al- 
lows us to cither write F = Vf(n,T) or F = Nf(v,T), 
with v = V/N — n _1 , such that / = nf. Hence 
both forms carry the same physical information. From 
the Helmholtz free energy we can obtain the chemical 
potential as p(n,T) = (df/dn) T and the pressure as 

p(n,T) = —(j)f/dvj . Therefore, p = p(n,T) and 

p = p(v,T) carry the same information. The laws of 
thermodynamics require that both the chemical poten- 
tial and the pressure are single valued, that is, for any 
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given values of n and T, or v and T, there must only 
exist a single value of p or p, respectively [24j . 

Therefore, a first thermodynamic condition that must 
be satisfied during the cosmological Bosc-Einstein Con- 
densation process is the continuity of the pressure at the 
transition point. With the use of Eqs. (fT4|) and (|25|) the 
continuity of the pressure uniquely fixes the critical tran- 
sition density p x from the normal dark matter state to 
the Bose-Einstein condensed state as 



u 



c 2 a 2 m x 
2nh 2 L ' 



(31) 



The numerical value of the transition density depends 
on three unknown parameters, the dark matter particle 
mass, the scattering length, and the dark matter particles 
velocity dispersion, respectively. By assuming a typical 
mass of the dark matter particle of the order of 1 eV (1 
eV = 1.78 x 10~ 33 g), a typical scattering length of the 
order of 10~ 10 cm, and a mean velocity square of the 
order of (if 2 ) = 81 x 10 14 cm 2 /s 2 , the critical transition 
density can be written as 



x 10 



-21 



10 



-33 



3 x 10- 6 

la 



10 



-10 



g/cm 3 . (32) 



The critical temperature at the moment of Bosc- 
Einstein condensate transition is given by [H-El 



T 



2nh 2 



C(3/2) 2 /3 m V 3 fcs V* ; 
(27r?i 2 ) 1/3 c 4 / 3 (a 2 f 3 m x /3 



C(3/2) 2 / 3 fci 



/ 



2/3 



(33) 



where £(3/2) is the Riemann zeta function, or 

1/3 



6.57 x 10 3 x 



10 -33 



3 x 10" 



2/3 



la 



10" 10 cm 



-2/3 



K. (34) 



The critical pressure of the dark matter fluid at the 
condensation moment can be obtained as 



P cr = 1.04 X 10" 



<7 



3 x lO- 6 

3 / , \ -1 



X 



In 



10" 10 cm 



dyne/cm 2 . (35) 



The critical value a cr of the scale factor of the universe 
at the moment of the condensation can be obtained from 
Eq. (J3TJ) as 



<W«0 - (-^2-J 



%2i o \ 1 /3(l+o' 2 ) 



c 2 <r 2 m 3 



(36) 



giving for the critical redshift of the transition the ex- 
pression 

/o *2, o N -i/3(i+<t 2 ) 

1+Z cr = 2 2 3 ■ ( 37 

By using the adopted numerical values of the constants 
we obtain for the critical scale factor and for the critical 
redshift the values 



a-cr/ao = 8.17 x 10 4 x 



10^33 g 
2 \ -l/3(l+a 2 ) 



3 x lO- 6 

, V l/3(l + cr 2 ) 

vn. \ 



10- 10 cm 



(38) 



and 



1 + z cr = 1.22 x 10 3 x 



10 -33 g 
2 x l/3(l+<r 2 ) 



(1+. 2 ) 



X 



3 x 10- 6 

/ \ -l/3(l+<r 2 ) 

1Q- 10 cm 



(39) 



respectively. 



B. Cosmological evolution during the 
Bose-Einstein Condensation phase 

During a first order phase transition, the temperature 
and the pressure are constants, T — T cr and P = P cr , 
respectively. The entropy S = sa 3 and the enthalpy 
W — (p + p) a 3 are also conserved quantities. After the 
beginning of the phase transition the density of the dark 
matter p x (t) decreases from p x r (T cr ) = p x or (when all 
the dark matter is in a normal, non-condensed form) to 
Px (T cr ) = P X EC '> corresponding to the full conversion of 
dark matter into a condensed state. It is convenient to 
replace p x (t) by h(t), the volume fraction of matter in 
the Bose-Einstein condensed phase, which is defined as 



h(t) 



Px (*) ~ PT 

n BEC n nor ' 



(40) 



Therefore the evolution of the dark matter energy density 
during the transition process is given by 

Px (*) = Px EC h(t) + P7 r [1 - Ht)] = Px° r I 1 + n xKt)\ . 

(41) 

(42) 



where we have denoted 



n BEC _ nor 
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At the beginning of the Bose-Einstein Condensation pro- 
cess h(t cr ) — 0, where t cr is the time corresponding to the 
beginning of the phase transition, and p x (t cr ) = p x or . At 
the end of the condensation h (t b ec) — T where tsEC is 
the time at which the phase transition ends, correspond- 
ing to p x (tBEc) = P X EC ■ F° r t > tBEC the universe 
enters in the Bose-Einstein condensed dark matter phase. 
From Eq. © we obtain 



1 



(p eec - pT) h 



1 rh 



a 3 p™ + P cr / C 2 + {pBEC _ p nor} h 3 1 + r/l ' 

(43) 

where we have denoted 



n BEC _ nor 



P™ r + Pcr/C 2 1 + P cr /p™ r C 2 



(44) 



where we have used the initial condition h (t cr ) = and 
we have denoted a cr = a (t cr ). At the end of the phase 
transition the scale factor of the universe has the value 



clbec = a{t B Ec) = a cr {l +r) 1/3 . 



(46) 



The Bose-Einstein Condensation modifies the overall ex- 
pansion rate of the universe. Consequently, during the 
phase transition, the pressureless baryonic matter and 
the radiation density evolve according to 



Pb 



Pb.a 



(a cr /a a y 



[l + rh(t)],t€(t cr ,t B Ec), (47) 



and 



Since p% EC < p™ r , generally r < 0, r G (-1,0), and 
n x < 0, respectively. Eq. (14*31 immediately leads to the 
expression of the scale factor of the universe during the 
Bose-Einstein Condensation phase as 

a(t) =a cr [l + rh{t)]- 1,Z ,t G {t cr ,t BE c) , (45) 



Prad = f fr ° J '°, [1 + rh(t)f /3 , t G (t cr , t B E C ) , (48) 

(acr/ao) 



respectively. The evolution of the fraction of the con- 
densed matter h(t) during the Bose-Einstein Condensa- 
tion process is described by the equation 



¥=-*(-+» 

(XT \r 



n 



I, A) 



(a cr /a Q )' 



(l + rh) + 



fi 



rad.O 



(a cr /a )' 



(1 + rh) 4/3 + fi x , wr (1 + n x h) + fi A , 



(49) 



where we have introduced a dimensionless time variable 



Hot, and we have denoted fL 



7po 



r/ P ™ r C 2 



a « 1, from Eq. 



, re- 
fill 



spectively. Since P c , 
it follows that r w n x , an approximation we adopt in 
the followings. If the energy density contribution of the 
radiation to the total energy density of the universe can 
be neglected, Eq. (j49j) can be integrated exactly to give 



h(t) = 



— X 

rlttr 

i + fi c 



id 

x exp [— (3fi A /in) (t — t cr )) 



1 - Slcond x exp [ 

rfltr 



(3fl A /t H ) (t - t cr )] 



(50) 



where we have denoted 



fi, 



n 



6,0 



(a C r/a y 



fix 



and 



fi 



cond 



y/Qtr + ^A — 



fi. 



fi A ' 



(51) 



(52) 



respectively. 

The time interval necessary to convert the entire exist- 
ing normal dark matter to the Bose-Einstein condensed 
phased is given by 



At r 



tBEC — t cr 



In- 



[\/fitr ~ 


hfi A " 


-fi A ] 


•y/fitr 


(1 + r)- 


f fi A + fi A 


[\/fitr " 


hfi A " 


^fi A ] 


V fitr 


(1 + r)- 


■f fiA — fi A 



X t H . 



(53) 



If the dark energy density can also be neglected, the volume fraction of matter h(t) in the Bose-Einstein con- 
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densed phase evolves as 



h(t) 



1 -2 



1 + -ifoV"tr — *c 



(54) 



In this approximation the conversion of the entire dark 
matter of the universe to the BEC phase is given by 



(1 + r) 



-1/2 



1 



H 



'Cttr 



(55) 



By assuming the standard values l a = 10~ 10 cm, m = 1 
eV and a 2 = 3 x 10~ 6 we obtain Q tr = 5.02 x 10 8 , giving 
for the total condensation time an approximate value of 



At 



cond 



= 1.31 x 



(1 + r) 



-1/2 



- 1 



x 10 



13 



(56) 



The time evolution of h(t), as given by Eq. (14U)) is 
represented, for several values of r, in Fig. [TJ 



1.0 



o.s 



0.0 



0.4 



0.2 



0.0 







/ / / 




1 t / 




/ y 








/ '' ' ' 




/ / 








/ / / ' ' 




: / / , / y 




/ / / / / 




/ / y 










: W' 




- >y 









1.x 10 



2.x 10 



3.x 10 
t/t H 



4.x 10 



5.x 10 



FIG. 1: Time evolution of the condensed dark matter fraction 
h(t) for different values of r: r = —0.10 (solid curve), r — 
—0.15 (dotted curve), r = —0.20 (dashed curve), r — —0.25 
(long-dashed curve), and r — —0.3 (ultra-long dashed curve), 
respectively. 



C. The Post-Condensation phase 

At t = tsEC and clbec = a cr (l + r) -1 / 3 , all the dark 
matter in the universe has been converted to a Bose- 
Einstein condensate phase. At this moment the cosmo- 
logical density of the dark matter is given by 



BEC 



C - ^%ii±^ , (57) 

uo (a cr /a ) -po x (l + r) 

I 



where the constant po x can be represented as 



Po x 



1.63 x 10~ 15 x (/ a /10- 10 cm) (to/IO^ 33 g) 



-3 



1 + 1.63 x 10- 15 x (/ a /10- 10 cm) (m/10" 33 g) 3 ' 

(58) 

The condition of the positivity of the density imposes 
the constraint (a cr /ao) > (1 + r) 1 ^ 3 on the model 
parameters. 

The equation determining the time evolution of the 
scale factor of the dark matter in the Bose-Einstein con- 
densate phase is given by 



d(a/a ) 
dt 



a/a 



yj(a/a f - po x 



(59) 



where we have denoted 
Use = 



^v.O 



1 + p x ,auo/c 2 '' 
and can be integrated exactly to give 



(60) 



x/QbeHo (t — C) 




arctan 



, (61) 



where C is an arbitrary constant of integration, which can 
be determined from the condition a = as EC &tt = tsEC- 
Thus we obtain 



C — tsEC — 



Pox 



'B£n 



2 / / a-BEC 

3VTIb~e~H V V «o 

-77 arctan 



Pa x - 3 x 



{a-BEc/ao) — P0x 

POx 



(62) 



In the case of a universe filled with dark energy, radia- 
tion, baryonic matter with negligible pressure, and Bose- 
Einstein condensed dark matter, respectively, the time 
evolution of the scale factor is given by the differential 
equation 



1 d(a/ao) 
(a/ao) dt 



a 



b,() 



n 



BE 



(a/a ) 3 (a/a ) 4 (a/a ) 3 



□a, t > tBECi 



POx 



(63) 
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which must be integrated with the initial condition 
a(tBEc) — fl (0) = as ec- The time evolutions of the 
scale factors for universes containing BEC dark matter 
are represented, for different values of the BEC parame- 
ter p 0x , in Fig.H 




-2.5 ; : 
-3.0 ; 1 

0.0 0.2 0.4 0.6 0.8 1.0 

FIG. 2: Time evolution (in a logarithmic scale) of the scale 
factor of the universe in the Post-Bose-Einstein Condensation 
phase, for cibec = o(0) = 8.1 x 10 -4 (z « 1200), and for 
different values of po x : po x = 10~ n (solid curve), po x = 
5 x 10 -11 (dotted curve), po x = 10~ 10 (dashed curve) and 
po x = 5 x 10~ 10 (long dashed curve), respectively. 

The presence of the condensate dark matter changes 
the global cosmological dynamics of the universe in 
the post-condensation phase, and the magnitude of the 
changes increases with the increase of the BEC parame- 
ter p 0x . 

IV. DISCUSSIONS AND FINAL REMARKS 

The Bose-Einstein Condensation process has been ex- 
tensively studied (especially for the case of scalar fields) 
in the physical literature. However, many important 
questions have yet to be answered. For example, one 
would like to know more exactly if the condensation pro- 
cess of the dark matter takes place (almost) instanta- 
neously, or it evolves during a long time interval. Many 
details of the thermodynamics of the Bose-Einstein Con- 
densation are not yet conclusively understood. Even the 
order of the phase transition (crossover?) is still a matter 
of debate. And, of course, from a cosmological point of 
view, one would like to know more exactly in what cos- 
mological epoch the condensation did occur, what where 
the values of the cosmological parameters at that time, 
how long the transition lasted, and what were its im- 
plications on the global evolution of the universe. In the 
present paper we have tried to give some preliminary and 
qualitative answers to these questions. 

Once the Bose-Einstein Condensation occurs, dark 
matter becomes a mixture of two phases, the normal and 
the condensed phase, respectively. From a general ther- 
modynamic point of view the coexistence of these two 



phases requires the continuity of the pressure of the two 
phases at the beginning of the condensation (phase tran- 
sition). This thermodynamical condition uniquely fixes 
the value of the dark matter density at the condensation 
moment, as well as the values of all the other thermo- 
dynamical parameters (temperature and pressure). The 
concrete numerical values of the condensation parame- 
ters depend on the mean square velocity of the normal 
dark matter a 2 , and of the scattering length l a and mass 
m x of the dark matter particle. Since the values of a 2 , 
l a and m x are not known for the dark matter particles, 
or are very uncertain (like m x and cr, for example), it is 
difficult to predict the exact cosmological moment of the 
Bose-Einstein condensation, and the corresponding val- 
ues of the cosmological parameters. However, by adopt- 
ing some "standard" numerical values, one can obtain 
a qualitative picture of the transition. Thus, the mass 
of the dark matter particle was assumed to be of the 
order of 1 eV [7], l27| > while the mean velocity of the non- 
relativistic dark matter particles was taken to be of the 
order of 900 km/s, which is somewhat larger than the 
values inferred from the study of the galactic dark mat- 
ter halos [13, The numerical value of the scatter- 
ing length l a for dark matter particles is very uncertain. 
As can be inferred from present experiments in ultracold 
gases, the scattering length is a quantity that determines 
the thermodynamic state of the gas. Terrestrial experi- 
ments on Bose-Einstein atomic condensates give a value 
of 10~ 7 cm [IH1]. For simplicity we have adopted for this 
parameter a value of l a = 10~ 10 cm. 

The general analysis of the condensation process shows 
a distinct three-phase cosmological history of the uni- 
verse. In the first phase, the universe evolved according 
to the standard A Cold-Dark Matter (ACDM) model. 
The condensation process of the dark matter did start 
when the temperature of the bosonic gas did fall below 
the critical value, and the equality of the pressure al- 
lowed the two phases to coexist simultaneously. With 
the adopted numerical values of the dark matter pa- 
rameters the condensation began at a redshift of around 
z = 1200. The condensation process took place grad- 
ually, the normal and condensed phases coexisting for 
around 10 6 years. To describe this transition period we 
have introduced the fraction h(t) of the condensed dark 
matter, and studied its time evolution. After the tran- 
sition phase the universe entered in the condensed dark 
matter phase. Depending on the moment of its occur- 
rence in the cosmological history, the presence of the 
condensate dark matter and of the Bose-Einstein phase 
transition could have modified drastically the cosmolog- 
ical evolution of the early universe, as well as the large 
scale structure formation process. Of course, even small 
changes in the numerical values of the parameter set 
(cr 2 ,Z a ,m x ) could lead to significant changes in the val- 
ues of the critical redshift or duration of the condensation 
phase. 

Immediately after the phase transition the presence of 
a Bose-Einstein condensate may drastically change the 
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cosmological dynamics of the universe. As one can see 
from Eq. if for some time interval (a/a ) is very 

close to pj^ 3 , the cosmological dynamic of the universe 
is determined by the condensed dark matter. In this 
case the condensed dark matter energy density p BEC — 

Qbe/ (a/oo) 3 — Pox i s very large, and dominates all the 
other cosmological energy terms. When this condition is 
fulfilled, the expansion of the universe enters in an accel- 
erated phase. For the time interval for which p BE can be 
approximated as a constant, p BE rj constant, the expan- 
sion of the universe can be described as a de Sitter one. 
The presence of an accelerating expansion period after 
the phase transition can also be seen from the analysis of 
the deceleration parameter q = —ad/ a 2 = d{\/H)/dt — 1 

1/3 

of the universe. For (a/ do) — > p 0x , the scale factor of 
the universe with Bose-Einstein condensate dark matter, 
given by Eq. (|62p , can be approximated as 



a 
f'o 



Cf + Pox 



1/3 



(64) 



During this phase the deceleration parameter q is given 
by 



1 

1 ~ n 



2p, 



■ox 



2 zn BE H%(t-cy 



(65) 



During the time interval given by t < C + 
a/ (4/3)po x /£IbeHq the deceleration parameter satisfies 
the condition q < 0. However, in the limit of large times, 
the expansion of the universe decelerates, with the decel- 
eration parameter given by q ~ 1/2. 

Therefore, the Bose-Einstein condensation of the dark 
matter may also provide the explanation for the observed 
recent acceleration of the universe, which could be just 
a transient phenomenon related to the recent entering 
of the universe in the Post-Bose-Einstein condensation 
phase. Hence BEC could also provide an alternative, 
physically realistic solution to the dark energy problem 
(for reviews of the dark energy problem see [8]). However, 
such a model would require that the cosmological BEC 
phase transition took place very recently, at z « 1 — 2, on 
a global cosmic scale. Presently there is no astrophysical 
evidence for such a scenario. 

A central problem in the theory of the Bose-Einstein 
condensation is the nature of the dark matter particle, 
as well as its physical properties. The core/cusp prob- 
lem can be solved if the dark matter is composed of ul- 
tralight scalar particles with mass m = 10~ 22 eV, ini- 
tially in a (cold) Bose-Einstein condensate [11]]. The 
wave properties of the dark matter stabilize gravitational 
collapse, providing halo cores and sharply suppressing 



small-scale linear power. Ultra-light dark matter parti- 
cles with masses of the order of m ~ 10 -23 eV and with 
Compton wavelengths of the order of galactic scales have 
temperatures of about 0.9 K, a temperature substantially 
lower than the temperature of CMB neutrinos. There- 
fore such a mass would guarantee that Big Bang Nu- 
cleosynthesis remains unaffected. In addition, the tem- 
perature is consistent with WMAP observations without 
fine-tuning [28j . It was shown that axions with mass 
to « 6 x 10" 6 eV 10 12 GeV// also form a BEC, where 
/ > 10 9 GeV This would imply an axion mass of 

the order of m f=s 10~ 3 eV. The axions were postulated 
shortly after the standard model of elementary particles 
was established to explain why the strong interactions 
conserve the discrete symmetries P and CP (for more 
details see [l4[ and references therein). A cosmological 
bound on the mass of the condensate particle can be ob- 
tained as to < 2.696 (g d /g) (T d /T cr ) 3 eV 0, where g is 
the number of internal degrees of freedom of the particle 
before decoupling, g d is the number of internal degrees 
of freedom of the particle at the decoupling, and Td is 
the decoupling temperature. In the Bose condensed case 
Td/T c < 1, and it follows that the BEC particle should be 
light, unless it decouples very early on, at high temper- 
ature and with a large g^. Therefore, depending on the 
relation between the critical and the decoupling temper- 
atures, in order for a BEC light relic to act as cold dark 
matter, the decoupling scale must be higher than the 
electroweak scale [27|. On the other hand, from Eq. ((T|) 
for l a 1 fm and R ss 10 kpc, it follows that the typical 
mass of the condensate particle is of the order of to ~ 14 
meV. For l a « 10 6 fm, corresponding to the values of a 
observed in terrestrial laboratory experiments, m w 1.44 
eV. These values are consistent with the limit m < 1.87 
eV obtained for the mass of the condensate particle from 
cosmological considerations @ . 

A better understanding of the numerical values of 
the Bose-Einstein condensation parameters (scattering 
length and dark matter particle mass), would be very 
helpful in obtaining accurate cosmological conclusions in 
the framework of the BEC model. Such an advance may 
also provide a powerful method for observationally test- 
ing on a cosmological scale the theoretical predictions of 
the Bose-Einstein Condensation model, and the possible 
existence of the condensed dark matter. 
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